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Abstract: We determine coherent states peaked at classical space-time 
of the Schwarzschild black hole in the frame- work of canonical quantisation of 
general relativity. The information about the horizon is naturally encoded in 
the phase space variables, and the perturbative quantum fluctuations around 
the classical geometry depend on the distance from the horizon. For small 
black holes, space near the vicinity of the singularity appears discrete with 
the singular point excluded from the spectrum. 

1 Introduction 

Quantum Black holes are still extremely mysterious, even after a decade of 
optimism, which arose with the microscopic counting of Entropy. What is 
the horizon in the microscopic picture? What happens to the singularity 
at the centre of the black hole? What is exactly Hawking radiation, and 
where is the end point, in a pure or mixed state? There is an immense 
amount of literature ([1, 2, 3, 4]) trying to address these issues, using various 
approaches, and reviewing them all is impossible here. But what emerges 
from them, is that we need a microscopic theory, which is equipped enough 
to address non-perturbative gravity, i.e. not a theory of gravitons, but a 
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theory incorporating all the non-linearity associated with strongly gravitating 
systems like the black hole, also the quantum states must have a suitable 
classical limit so that one can identify continuum geometry. In otherwords, 
a non-perturbative quantum theory, where one can define h — > 0, without 
making the gravitational fields weak. String theory as it stands today, does 
not have a non-perturbative, back ground independent description so far. 
Other approaches to quantisation, involving Path-integrals and discretization 
[5], may have useful models, but are yet unexplored. This brings us to 
a very interesting construction, in the framework of canonical quantisation 
of general relativity, of coherent states, which supposedly in any quantum 
theory encode the information of classical objects. Thus, we have a theory of 
non-perturbative quantum gravity, which is back ground independent, and a 
construction which allows one to look for the classical limit. Hence one should 
immediately apply them to black holes and investigate the appearance of the 
horizon and singularity in their framework. The task is thus two prong, 
trying to understand what coherent states are in the context of gravity, and 
their application to extract information about a classical black hole. 

1.1 Coherent States 

In quantum mechanics, coherent states appear as eigenstates of the an- 
nihilation operator, and in these the minimum uncertainty principle, i.e., 
AxAp = h/2 is maintained by the position and momentum operators. This 
definition has a generalisation in terms of a 'coherent state transform' defined 
thus: Given a space L 2 (R n ), then there exists a map or an integral transform 
to the space Ti(C n ) of holomorphic functions on C n [6]. In other words: 



Where A(z,x) = ( 27r )n/4 ex P [~\^(^ z i + x \ ) + ^ x i z i] , i n case of i — 1, 
this is same as the simple harmonic oscillator coherent state wave function. 
The above also has a more convenient representation 



With p t being the heat kernel on R n , and pt(z) is the analytic continuation 
of p t to C n and t denotes the width of the distribution. This is also called the 
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Segal-Bargmann transform. This definition was generalised for an arbitrary- 
gauge group in [6], with the following: 



A,( 9 , x) = ^ (3) 
VPt( x ) 

with x belonging to a compact gauge group K, g belonging to a complexifica- 
tion of the gauge group, G, and the pt{x) being the heat kernel on K defined 
with respect to an appropriate measure. Once the A t (g,x) are determined, 
they will serve as coherent states in the 'position representation' with width 
y/t, and there will be a corresponding one in the 'conjugate momentum' rep- 
resentation. It is this definition, for the coherent states which can be applied 
to gravity, as it appears as a SU(2) gauge theory in the canonical framework. 
For further work on Coherent States see [7] . 



1.2 Coherent States for Gravity 

In the canonical framework, space-time is separated as R x E, where E is a 
3-dimensional manifold, with a given topology. The two physical quantities 
relevant for this spatial slicing are the triads which determine the induced 
metric and hence the intrinsic geometry of E, and the extrinsic curvature 
which determines how the slices are embedded in the given space-time. The 
triads have an internal degree of freedom associated with the tangent plane. 
This 5*0(3) symmetry on the tangent plane, is isomorphic to SU(2), and 
one can combine the variables to give a connection which transforms in the 
adjoint of the SU(2) gauge group. The precise definitions are: 

Ai^ri + PK^e" 1 , (3E?E bI = detqq ab (4) 

where a denotes the spatial indices 1, 2, 3 and I the internal SU(2) index 
also taking values 1,2,3. q a b is the induced metric on the three manifold, 
the e b j are the corresponding triads, is the spin connection associated with 
them. K a i, is the pull back of the extrinsic curvature on to the spatial slice. 
The new variables, the connection A\ and the one form constitute what 
is known as the 'Ashtekar-Barbero-Immirzi' variables, with a one parameter 
ambiguity present in the definition called the Immirzi parameter (5. Details of 
this can be found in the reviews [8, 9, 10] and books [11]. The Einstein-Hilbert 
action written down in the above variables includes additional constraints 
associated with the lapse and the shift, as well as the Gauss constraints 
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associated with any Yang-Mill's theory. On the constrained space, the 'Phase 
Space' variables satisfy appropriate Poisson brackets. Since the quantisation 
has to be background independent, and non-perturbative, one takes the help 
of 'Wilson Loops' which are path-ordered exponentials of the connection, 
along one-dimensional curves. This leads to the introduction of abstract 
graphs F(v, e, J) which are comprised of the disjoint sets of vertices v, edges e, 
and a incidence relation / between them. e.g. A graph with every element of 
e is incident with two elements of v is a simple planar graph. The restriction 
on the graphs used in this formalism is the requirement of piecewise analytic 
edges. For details, see [8, 9, 10]. The holonomies along these edges constitute 
the configuration space variables. The quantum Hilbert space constitutes 
of Cylindrical functions, f(h ei (A), h e2 (A), h e3 (A), h ei (A), ..) where h ei (A) 
is the holonomy of the connection along the i th edge. This brings us to 
the coherent state transform within the framework of loop quantum gravity 
[14, 15, 16, 17, 18, 19], For each edge e, the holonomy h e (A) C SU(2), and 
hence a complexification of K=SU(2), would lead to the complex group G 
defined by Hall. The transform would take a function cylindrical over SU(2), 
to a function of the complexified holonomy. The transform will have the 
kernel, using the generalisation by Hall. There is a subtlety in the definition 
pointed out by B. Hall, where the actual transform used in gravity is the C t 
transform defined in Equation C, of Appendix of ([6]). 

C t {g,h)=Y,d«e- tK X«[9h- 1 ] (5) 

Where in the above, h is the element corresponding to the holonomy, g, 
the complexified group element, and the RHS is an expansion of the Heat 
Kernel on SU(2) based on Peter- Weyl decomposition. Given an irreducible 
representation n of SU(2) or its complexfication SL(2,C), d n is its dimension, 
A,,- the eigenvalue of the Laplacian, and Xn its character. Thus, the coherent 
state on each edge of the graph by above is: 

ri(g,h) = C t (g,h) (6) 

The coherent state for the entire manifold is then 

V t r (g,h) = l[ri(g,h) (7) 

e 

Thus given the above, it remains to find the complexified group element g 
for gravity. Usually the conjugate variable to A^, the densitised triad is used 
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to complete the phase-space, but for requirements of a anomaly free algebra, 
and gauge covariance, Thiemann has defined a momentum variable 



P 1 = Tr 

a N 



T% (j h p (p)E a h; 1 (p)e abc dS bc: ] K 1 



(8) 



The 'conjugate momentum' is associated with each edge 'e' and evaluated 
on a D — 1 'dual' surface, which comprise a Polyhedronal decomposition of 
the sphere. It does not depend exclusively on the densitised triad, but also on 
the connection through the holonomies h e , and h p as defined above. One thus 
needs a dual polyhedronal decomposition of S, made up of dim D-l surfaces, 
intersected by the edges transversely. The holonomies h e are holonomies 
evaluated on the edge from the starting point of the edge to the point of 
intersection. The h p are holonomies along edges confined to the dual surface, 
extending from a given point p to point of intersection. (T 1 is the generator) 
Now, this specially smeared operator, has an appropriate Poisson algebra 
with the holonomy variable h e (A). 



{h e , h e ,} = {h e , P/,} = tb^T'K {P/, P e J ,} = te IJK P e K 6 ee , (9) 

(t = n/a N ). This algebra can be lifted to commutator brackets, where now 
t = Zp/ajv and has a suitable representation. Using this and a suitable com- 
plexification scheme, Thiemann determined the complexified group element 
to be 

g = ^KiA) (10) 

Note the presence of an undetermined constant which was introduced 
to make the quantity P dimensionless. This quantity as noted in [20] has an 
analogy with the Harmonic oscillator frequency, and is a input in the theory. 

Thus, given a graph T, its piecewise analytic edges, the corresponding dual 
graph, the h e , P/ constitute the graph degrees of freedom. The coherent state 
is obtained as a function of these parameterized as h e = T , P = P / T/. 
The knowledge of what the classical values for these are once a graph is em- 
bedded in a space-time which is a solution to Einstein's Equation, the point g 
at which the state must be peaked is determined through Equation (10). As 
defined in Equation (5), the coherent state is completely known. Using the 
above definition for the coherent state, and it's gauge invariant counterpart 
obtained by group averaging method, in a series of papers, these coherent 
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states were demonstrated as having the following: They were peaked at ap- 
propriate classical values, they satisfied the Ehrenfest theorem, (expectation 
values of operators were close to their classical values, and one could obtain 
a perturbation series in t); the states are overcomplete with respect to an 
appropriate measure [15, 16, 17, 18]. Finally to make the definition complete 
(7), one must take Infinite tensor product limit to incorporate asymptotically 
flat space-time (also the case relevant here), where the Graphs are infinite. 
This was addressed in [18, 19]. So, modulo the dependence on the Graph, and 
the choice of corresponding dual surfaces, one can believe that one has found 
the coherent states for non-perturbative gravity. For a comprehensive dis- 
cussion, and other attempts to relate to semiclassical physics see [21, 22, 23]. 

1.3 For a Black Hole 

What would thus a coherent state for a given classical space-time mean? The 
states as mentioned earlier are defined on edges of a given abstract graph. 
This graph can also be embedded in the classical spatial slices, and the 
distances along the graph measured by the classical metric. The information 
about the classical holonomy along the embedded edge is encoded in the 
complexified group element g. The behavior of ip t as a function of h in the 
connection representation gives the quantum fluctuations around the given 
classical value g. The construction of the coherent state for a black hole 
would thus involve, choice of a suitable slicing of the black hole, embedding 
the chosen graph in the slice, and then evaluation of the classical holonomy 
and the corresponding momentum (8) along the edges. 

In this article we use the above outlined framework to determine the 
coherent states for a Schwarzschild black hole. The formalism of course can be 
extended to other black holes with multiple horizons, non-trivial maximally 
symmetric space-times like de-sitter space which has a cosmological horizon. 
Using a suitable choice of coordinates for the static black hole space-time, 
we obtain flat spatial slices, and embed spherically symmetric graphs on the 
slices. All these exersices are done to extract maximal information in the 
simplest possible settings, as apriori the task looks quite a difficult one to 
handle. Once the classical holonomies and the momenta have been evaluated 
along the embedded edges, one can write down the wavefunction for each 
edge. Since the position of the embedded edge can be measured by the 
classical metric, one can comment on how the quantum fluctuations behave 
as a function of distance from the center of the black hole. We find that 
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the wavefunction depends on the position of the horizon, and for a particular 
orientation of the radial edge, becomes a gaussian. This radial edge lies along 
the 9$ = n/2 plane , and we call it the 'ray of least resistance'. We follow this 
edge deep into the black hole using finer and finer dual surfaces, and calculate 
the corrections to the geometry. The idea in this article is to evaluate the 
corrections to the momentum and holonomy operator, and since most other 
quantities can be defined in their terms, the quantum fluctuations to other 
operators can be determined using these. The results we find are interesting, 
but not dramatic. The momentum operator for the radial edge has the 
information about the horizon, as expected as the dual surfaces are r = const 
and hence lie transverse to a freely falling particle. Depending on the size of 
the edge, oscillations set in as soon as one crosses the horizon, with decreasing 
amplitudes. The tidal forces acting on the spherical surface distort it, and 
oscillations arise due to that. The frequency increases, but the amplitude 
of oscillations are as small as 1CT 9 . The P are measured with respect to a 

given length scale a N and this is what determines the 'semiclassicality of the 

i 2 

wavefunction'. The semiclassicality parameter is t — and hence the bigger 
is with respect to the Planck Length, the more is the wavefunction peaked 
at the classical value. However as we discuss later on the value of the classical 

r 2 

momentum falls sharply once inside the horizon, and is proportional to 
Hence to keep the amplitude of oscillations regulated as one falls inside the 
black hole to 'semi-classical values', it has to take values much smaller than 
the black hole size. Optimizing leads to an interesting observation that 
for black holes which are slightly bigger than Planck size, the semiclassicality 
parameter is around 1, preventing a perturbative expansion around a given 
classical value. 

For larger black holes, everything is smooth across the horizon, as one 
is in a freely falling reference frame. Note at this stage, one is calling the 
position of the horizon, as the point at which r = r g , i.e. the classical 
value. This brings us to the derivation of the corrections to geometry, in the 
present situation, they look perfectly in control. We calculate the corrections 
in the P 1 component, for the radial edge, which lies along the path of least 
resistance. For large P, the corrections go as t(P + l + 1/P), and for small P 
as t(P 2 + P + 1). Since P is well behaved, the corrections are not unbounded 
anywhere in the configuration representation. 

To examine the singularity one has to look into the details of the em- 
bedded graph, which taking into account the singularity of the connection, 
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cannot be embedded classically at r=0, but starts some distance away. The 
dual graph, however includes the r = point, as the densitised triad by itself 
is not singular at r = and there is no harm in including the same. Thus 
despite the fact that one starts in configuration space without including the 
singularity, the dual graph when embedded has the knowledge about the sin- 
gular point. Thus, in the vicinity of the singularity it is correct to work in 
the Electric field representation, and that is what we resort to. 

After examining the semiclassical expansion around a given classical value, 
one can extend the argument to include values of h, P whose value is ~ t, 
(Note t is a dimensionless number). To do that, and also examine the space 
near the vicinity of the singularity, one goes to the electric field or momen- 
tum representation, and finds that the peak for small black holes, is restricted 
to take certain quantised value of the momenta. (Note that it is the dual 
graph, which can include the singular point, as it is does not diverge there) 
For large values of P, and large enough black holes, the classical continuum 
geometry is recovered. However, for smaller black holes, as one approaches 
the singularity, the momentum is peaked only at certain discrete values given 
by (j + l/2)t. The important point is the presence of a non-zero 'remnant' 
which prevents the momentum operator to be ever peaked at the singular 
point. So in conclusion, quantisation in terms of the above operators, should 
not see the singularity. One could question the legitimacy of such a strong 
claim given that one is working with the simplest graph, but this is related 
to a property of every single radial edge, which can be part of any graph. 

In the next section we begin by slicing the Schwarzschild space-time, then 
determining the appropriate classical variables. In the 3rd section, we work 
with the simplest possible graph for this spatial slicing, and embed it. After 
evaluating the classical group element, we investigate a single radial edge, 
which is reminiscent of the straight edge use in loop quantum cosmology 
[25, 26, 27]. There the justification used is the isotropy of space, and here, 
one uses the inherent spherical symmetry. In the next section, the Coherent 
state obtained using this edge is investigated in details. In the fifth section we 
include the angular edges of the graph, and briefly comment on the coherent 
state for the entire graph. Section six involves a conclusion and points for 
further discussion. 
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2 The Classical Preliminaries 



The coherent states as defined are on a spatial slice at a fixed time. So, here 
we take the static Schwarzschild metric, and look for a suitable slicing, such 
that the spatial slices are flat, and time has a definite direction. The usual 
Schwarzschild metric in static coordinates has the following form: 

ds 2 = - (l - — ) dt 2 + . dT \ , + rW (11) 
V r) (1-i) 

with time t, radial coordinate r. It is well known, that this coordinate 
chart covers only r > 2GM(= r g ). And would not be a very convenient one 
to work with, neither are the t = const slices flat. The other very popular 
coordinates are the Kruskal coordinates and then the Eddington-Finkelstein 
coordinates, both of which are suitable for studying null slices. The choice of 
slicing for addressing this problem should be compatible with the canonical 
frame work, where time is distinguished, and hopefully cover the entire Black 
hole space-time. One obvious choice is then the following metric, where 
space-time is foliated by a set of proper-time observers hurtling down through 
the black hole horizon into the singularity (This set of coordinates are also 
called the Lemaitre coordinates). 

The metric written in proper-time coordinates is: 



, 9 dR 
ds = + + 



2r 



\(R-r) 



4/3 

r 2 J 3 (d9 2 + sin 2 Odcf 2 ) (12) 



with r denoting the proper-time, and R the spatial direction. If one looks 
at the coordinate transformations, one finds the following relations [12], 



' dr=(dR-dr) (13) 



2GM 



( 2r \ 2/3 

dt = ir-f d R, f=( W ^ T) ) ("I 

From (13), it is clear that r = const surfaces are at 45 degrees in the R — t 
plane, and the coordinates are smooth across r = 2GM. Also, the singularity 
of the Black hole is at R = t. The constant r slices are flat, The slices have no 
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intrinsic curvature, and this becomes evident by making the transformation 
(at constant r = r c ) 

(h ' ' (15) 



dR 



1/3 : 



the metric is flat dr 2 + r 2 d£l 2 . The extrinsic curvatures given by Kij = 



K 



RR 



2r n 



-(R-r)\ 



5/3' 



2r, 



-(R-r) 



1/3 



r g , K(p(p 



K BB sin 2 9. 



(16) 

As expected the extrinsic curvature is singular at the centre of the black 
hole, showing that the slicing becomes more and more curved in the interior 
of the horizon. However, there is apriori no reason to treat the horizon as a 
boundary as the metric is perfectly smooth across it. However one must add 
a note of caution to the above, as the coordinates at present are not complete, 
in particular, for geodesies starting very far in the past, they only graze the 
horizon, never crossing it [13]. For the article in hand, this is not a major 
problem, as we are not addressing global issues like the presence of 'event 
horizon'. The constant r slices across some finite range in time are enough 
for us to build the coherent states, and the information about the horizon is 
in the form of the apparent horizon, or the outermost trapped surface. The 
pullback of the above to the spatial slice is very simple to derive from the 
above: 

Krr = \ \ —^4- I . Khh - * fn\, K,„, = - ^[W g sin 2 9 (17) 



-3/2 



Thus, at a constant proper-time r c , the classical spatial slice is flat as stated 
earlier, and one proceeds to a determination of the coherent state wave- 
function on a given slice S Tc . Throughout our work, our only input is the 
classical metric, and we expect the formalism to have answers regarding the 
geometry and the quantum corrections. Thus, to use the formalism devel- 
oped in [15, 17], one needs a graph to be embedded in the given spatial-slice, 
and then evaluation of the h e , P* with the classical metric in the embedded 
graph. The coherent states will be peaked at these classical values, and it is 
the quantum corrections and the nature of the wavefunction as a function of 
h e , P/, which crucially depend on the classical metric, will give us the answer 
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to the question we are asking: Does the quantum geometry know about the 
horizon, the singularity, existence of trapped surfaces etc etc. Thus as a be- 
gining, one needs to determine what the canonical variables are for the above 
slicing. As given in Equation (4), one can determine the connection and the 
densitised triad by using Equation(17), and the fact that the induced metric 
is flat. There is a gauge freedom in choosing the triads, and one can choose 
them to be diagonal in the indices, thus 

e\ = \ e e 2 = - e \ = ^— (18) 

Using the above and (17), one gets 



Al = -(3K rr e rl = -^ , A 2 e = pJ^, 4 = /3 A Msin0 (19) 





The one forms or the electric fields are: 

r 2 sin 9 * r sin 6 + r . . 

E\ = El = — EZ = —= 20 

1 v 7 ^ 

Thus once the classical variables have been determined, one can proceed 
to the abstract graph space, which when embedded will give us the infor- 
mation about the holonomy and the dual momentum. Note that there is 
additional contribution to the gauge connection from the spin-connection 
from the spherically symmetric sector of the metric, which is independent of 

T 9- 



3 The Simplest Graph 

This brings us to the question of what graph we must employ here, so that 
one can embed it in the given spatial slice. Since the metric has spherical 
symmetry, one would presume the use of a spherically symmetric graph. This 
is motivated from the work on loop quantum cosmology where the reduced 
Hamiltonian constraint forces the choice of a straight edge [25, 26, 27]. What 
would a spherically symmetric graph look like? Edges which start out from 
some central point like spokes in a wheel, and vertices along fixed distances 
along the spoke. The vertices will constitute a cycle, be inter connected by 
angular edges such that the last edge ends on the first vertex. Given the fact 
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that apriori nothing is clear about what the coherent state is going to look like 
for the Schwarzschild black hole, one can start with this simple graph indeed. 
Enlarging the graph to include more and more edges is not a problem as one 
can consider one to be a subgraph of the other. Also, as defined in (7), the 
coherent state is obtained as a product of the coherent states for each edge. 
Hence introducing more edges contributes to the product. So as a begining, 
one starts with a planar graph, a extended version of the more well known 
'Wheel graph'. Which has n 3-valent vertices with the % th vertex adjacent 
to the i — j th vertex and the i + j th vertex in the list , and one special 
vertex also called the hub which is adjacent to all the other n vertices. This 
wheel graph can be extended to include more vertices which are adjacent to 
the previous n vertices constituting a cycle (The initial vertex is same as the 
terminal vertex). The ith vertex in this new cycle is adjacent to the % — j th 
vertex and % + j th vertex, and to the ith vertex of the previous cycle 2 . One 
can go on adding new cycles. Now, one can add an infinite number of cycles 
ending up with a countably infinite graph. However, since the embedding 
will involve the hub situated at the origin of the singularity, one excludes the 
hub from the configuration space graph. 




Figure 1: Spherical Graph 
2 For a introduction to Graph theory [24] 
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There are discussions on the exact nature of graph which shall sample 
maximum number of continuum points, [28], e.g, one can modify the graph 
in Figure(l). Here one has added some more vertices to the cycles, but the 
edges connecting them are no-longer be the 'spokes of a wheel'. This is as 
suggested in Figure (2). 




Figure 2: Non-radial edges 

However for our purposes we confine ourselves to the Wheel graph, as this 
is relevant to us. The embedding of these planar graphs in the 9 = c and the 
= c planes can give the entire sphere, just like a latitude-longitude grid. 
A polyhedronal decomposition of the sphere csn be obtained by determining 
the planar dual and evolving them along the <fi or 9 directions. This fixes 
the 'dual' surfaces uniquely. The other way of understanding this is that the 
planar dual is the intersection of the dual polyhedrons with the corresponding 
planes in which the planar graphs are embedded. The planar dual graph 
corresponding to the graph described earlier is again a set of adjacent cycles, 
and the hub however reappears here. 

Question is now, how does one embed this radial spoke into the above 
spatial-slice. One point to note is that the r = constant slices for the 
Schwarzschild radial coordinate, are actually R — r slices of the Lemaitre 
metric. When one is emedding the radial spoke, in the above spatial slice, 
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Figure 3: Dual Spherical Graph 



which is a r = constant slice, and measuring distances in dr, defined in (15), 
one has to keep r c constant, i.e. stay in the same slice. Once the graph is 
embedded, in the slice, one can determine the holonomy and the momentum 
(8). As per the definition (5), one chooses to define the complexified group 
element g on the embedding of the edge in the classical space time. This 
would involve the holonomy and evaluation of the momentum (8) on the 
dual surface. 

3.1 The Complexified group element 

The complexified group element g as stated in Equation (10) has information 
about the classical holonomy and the corresponding momentum P evaluated 
along the dual surfaces. We thus begin by evaluating the holonomy on the 
embedded radial edge. 

3.1.1 The Radial Edge 

By definition, the holonomy along a given curve (here being the radial edge) 
is the path ordered exponential of the connection integrated along the curve. 
In other words: 
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h e (A) = Pexp( J A^dx") (21) 

fb fb rb j n j n, 

/ dAi / d\ 2 ... / dA^— .... (22) 



Where A is the parameter along the curve. For a radial edge, embedded in 
the spatial slice, (A is identified with r), and due to the diagonal connection, 
the holonomy is: 



h r e (A) = 1 + ^—^ / A r dr + 



ri 

COS j T < — — -7- ^ ) — KT 1 Sill 






1/2 1/ I I 1/2 1/2 



Where r = f an d the SU(2) generator T* = — |<r\ The radial edge 
begins at the point r\ and ends on the point r 2 as measured by the metric 
on the sphere. Since only one component of the connection contributes, the 
subsequent terms commute, and one is left with the above summable series. 
Of the few things to observe in the above holonomy calculation is that due 
to the singularity of the integrand, one in principle must integrate from a 
region e distance away from r = 0. So the classical curvature singularity is 
manifest in the classical holonomy, and whether it is anyway different in the 
quantum expectation value is the thing to explore. Next, one must evaluate 
the corresponding conjugate momentum. This as defined in Equation(8), 
involves the choice of a dual surface, and then a choice of edges on the dual 
surface which one must integrate over to get the momentum. By definition, 
this dual surface has to be chosen such that the radial edge intersects the 
surface transversely at a point p, only once. At the point p, the orientation 
of the dual surface is in the direction of the edge e. All these conditions are 
satisfied by pieces of spherical surfaces centred at the origin of the radial line. 
The dual polyhedronal decomposition is thus made up of spheres, but now 
cutting the radial edges at the midpoint. This ensures that the faces of the 
Dual graph are intersected transversely by the edges, at precisely one point. 
The intersection of the dual 2-spheres with the plane of the graph gives the 
edges of the dual planar graph Figure(3). The Radial edges and their duals 
can thus be visualised as in the Figure (4). 
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Figure 4: Radial Edge and Dual Surface 



The h p (9, 0) are holonomies along a edge connecting a generic point p 
on this surface to the point of intersection. This choice of the edges on the 
dual sphere is again difficult to make on the spherical surface, without intr 
oducing any ambiguity (see Appendix 2). In the following, the edges are 9, 
lines, and the shortest distance along the graph edges as an attempt to fix the 
distances uniquely does not give a meaningful answer as given in Appendix 
2. In principle any generic curve connecting the point (9, 0) with the point 
of intersection p would also suffice. In this article, we confine ourselves to 
the = plane, and the dual surface, which has to be comprised of the 
2-dim surface is a tiny strip extending from = + e/2 to — — e/2 
(e being a very small number). The intersections of these dual surfaces are 
one-dimensional edges on the = plane, and these edges comprise the dual 
planar graph. A discussion on the full three dimensional momentum is given 
in Appendix 1, and this preliminary simplified 2-surface is sufficient to illus- 
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trate the emergence of the horizon in the variables, and then in the coherent 
state. The simplified 2-surface has edges along the curve 9 = aX. Where A 
is the parameter along the curve and a some undetermined constant. One 
can think of this as the equator, and the radial edge intersects this surface 
once at some point 6$. The holonomy from any generic point (9, 0) is eval- 
uated along this curve. Eventually one will integrate over the contributions 
to the holonomy. What we are actually doing in this article is suppressing 
the azimuthal angle 0, which as we show in Appendix 1 does not change the 
results drastically. The information that the surface is actually a two sphere 
is encoded in the connection, and the one form E, which is a densitised triad 
and includes ^/q, where q is the three metric. 
The holonomy on the very thin strip is thus: 



h(e,e ) = i+ / Agde + .... 



00 
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2^2 



(0o " 0) 
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[acr 2 + ex 3 ) sin 
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(25) 



(0 O - 0) 



with a = afly -f. Note now one has taken the contribution from the spin- 
connection of the spherically symmetric sector into account. We now station 
the surface at a radial coordinate r, and the edge intersects the surface at the 
halfway point. The length of the edge is 5. This gives the radial holonomy 
from Equation (24) with the initial point r\ = r — 5/2 and the final point 
r 2 = r + 5/2 for r 3> 5 to be h e (r) = cos (^) — ta 1 sin (^) . Using the above 
and (25) in (8), and one gets the following expression for the momentum 
variable : 



P; = —Tr 

a N 

Where 



a cos 
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(26) 



X(r) = [ h(9,<j)) *E h-\9,(f)) 
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(27) 
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(The gives a overall factor of e which is not shown) This when evaluated 
gives 



1 



X(r) = X^r)* 1 + X 3 (r) ^— V + X 3 J?— a 3 (28) 



with 
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(29) 
(30) 



Where a' = Now, one uses the above in (26), to get the following 

result for the momentum components 



P 1 



Mr) 
a N 



p2 

r 



pi 



a N 



, [asm (7a 3 ) + cos (7a 3 )] (31) 

v a 2 + 1 

3 [— a cos (7a 3 ) + sin (7a 3 )] 
V'n- • 1 



The most important thing to notice in the above is the appearance of the 
factor 1/(1 — a 2 ) in the above. This for (3 = l/(a) denotes the Schwarzschild 
factor, (a = 1 and for our purposes we set (5 — 1, it can be restored in 
the final results, and it's value fixed when comparing to known semiclassical 
quantities like the entropy etc) and when the edge is at r = r g , this factor di- 
verges, signifying a abnormality in the behavior of the function P 1 . However 
things are smooth when one takes the limit carefully due to the presence of 
sin(l — a!)Q' in the numerator, which also goes to zero at that point. (We 
subsequently replace the 9' = 9, and S/2r g = 7) Note, we did not start work- 
ing in a coordinate chart where the horizon was manifest, and it is a nice 
feature of the above construction that the horizon starts playing some spe- 
cial role. Before going into the discussion of what all these quantities mean, 
one important point to note is the 'graph' dependence of the quantities. The 
momentum crucially depends on the point at which the radial edge intersects 
the dual surface, and also the size of the edge, and the dual surfaces. Thus 
apart from the dependence on the classical metric, and hence the radial co- 
ordinate r at which the edge intersects the dual surface, the size is encoded 
in 5 and the width of the dual surface in 9. In some sense these are the graph 
degrees of freedom, and these additional parameters label the graph, and 
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its particulars. Before going into the construction of the coherent state, one 
must have a understanding of the quantity of P 1 for the given classical met- 
ric, and its behavior. To do that, we plot the following graphs, illustrating 
the individual components (in group space) and the gauge invariant quantity 
P = \/ Pi + P 2 + Pi- Also it is interesting to compare the quantity for the 
same graph embedded in flat space, with the above as one does not have a 
natural understanding of the momentum P e 7 . Thus the momentum variable 
(for the radial edge) will be 



^2 



PL ~ — (32) 

One sees a distinct difference between (32), and (32), as the first one has 
'non-abelian' nature, despite the same graph being embedded, and distinct 
new length scale in the form of r g . As one can see from (32), this is clearly 
something which goes as r 2 , whereas for the black hole, there is a modulating 
factor present which we investigate now. For the black hole, where what we 
call the P flat (r 2 ) dominates for < a < 1 or outside the black hole, and then 
around a = 1 the modulating factor P 1 jr 2 start playing a greater role, taking 
non-trivial oscillatory nature, with decreasing amplitude. Also, this behavior 
is not independent of the size of the graph, with the oscillations setting in 
later and later behind the horizon as the graph is made finer and finer. So, 
in some sense a = 1 is a transition region in these set of variables, and as we 
show in the later sections, it is the manifestation of the momentum variable 
P in the coherent wavefunction that makes the above analyses worthwhile. 
Let us start by examining the gauge invariant quantity P. 



P r = — X /Xf + Xl (33) 
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+ 2cos(2# ) 
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sin[(l - a')9] sin[(l + a')6} 
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The above is plotted in Figure (5) (In this and all subsequent Mathemat- 
ica generated figures, the angles 9, 9 appear as c, c', the variable a as a and 
the momenta as P's the main reason being Mathematica or mine inability 
to get greek symbols and superscripts in the exported .eps files). 

(One sets r 2 /aN = 200 as a test and 9 = tt/6 in the above.) What is 
evident from the above is that the fall off of the function as l/a 4 dominates 
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Figure 5: Gauge Invariant momentum 



for < a < 1, which is same as the behavior of the function as would be 
expected in flat space-time. Just to show, if one divided out the flat behavior, 
one would get a modulated graph as shown in Figure (6). Around a — 1, 
one has to take the limit carefully, and one gets the P r as a function of the 
angle at which the edge intersects the dual surface (6*0), as well as the size 
of the dual surface 9. Now, one must notice, that there is nothing drastic 
in the behavior at the point a = 1. This is as plotted as a function of the 
angle 9, 9 (c,c' in Figure (7)) at which the radial edge intersects the horizon 
in Figure (7). 

However, around a > 1 oscillations set in, depending on the size of the 
dual surface. If the dual surface is big enough (e.g 9 = ir/2) oscillations 
set in as near as a = 3, for finer and finer graphs, 9 = 7r/100, oscillations 
set in around a = 200. The reason for this is that in the argument of the 
sin[(l — a')Ti/N] functions, there are two competing factors 1 — a' and tt/N. 
As soon as 1 — a' > N, the oscillations set in. This is an indication of the 
inability of the surface to be spherical in the increasingly curved slice of the 
black hole. The geodesic deviations become larger and larger, creating the 
above situation. Below is the plot for 9 = 7r/100 Figure (8). 
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Figure 6: P r /Pfiat 

Similarly, one can obtain the plots of individual components as shown 
one by one below. 

This first figure shows the P r x component as a plot of < a < 10 denoted 
as the a axis and 9 denotes as the c axis. The value decreases progressively 
about a — 1, it has very low values indeed Figure(9. 1). 

This next graph shows the same momentum but with a now in the range 
(100,1000) Figure(9.2). The oscillations set in with decreasing amplitude. 
(Though towards the later range the graph appears flat due to lack of reso- 
lution) 

The plot Figure (10) is a graph of P r 2 set at 9q = 7r/6 and the size of the 
dual surface is as large as ir/2. The Oscillations set in as quickly as a = 3. 
A plot of P r 3 is also quite similar to the above. 
The few points to note in the above discussion are: 

1) The momentum P r is oscillatory inside the horizon. 

2) The amplitude is very small, around 10~ 9 with r 2 g ja^ = 200 , just inside 
the horizon and decreases progressively. Thus to prevent the momentum 
from reaching very low values, one needs to set = 10 _9 r 2 at least. The 
implications of this are discussed in detail in Section 4.2 

3) The finer the graph is later the oscillations start, indicating an association 
with the geodesic deviation forces which act on the embedded edges. 
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Figure 7: Momentum at r = r g 



3.1.2 Angular edge 

Now we come to the evaluation of the complexified group element for the 
angular edge. The angular holonomy is already evaluated along a 9 curve 
in Equation (25). Here we directly evaluate the angular edge momentum. 
Along the plane, it is again a small strip enclosing the radial line along the 
dual planar graph. 

We thus determine P$ , and this has the form: 



Where when one puts in the expressions for holonomies (24,25) , one gets 




-Tr[T% Y(6) h 1 } 



(34) 



Y(0) = Y 2 a 2 + Y 3 a 3 



(35) 



with 
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Figure 8: Finer Graph 





Figure 9: P\: Across the Horizon and Inside the horizon 



The above two integrals can be done, and yield 
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Figure 10: P 2 : Large size graph 
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In the above a + = 1/a 2 + 7 and a_ = 1/a 2 — 7. (7 = S/2r g ) and Ci and S'i 
denote cosine and sine integral functions respectively. The momentum is: 



Pe = — 

(In 



Y 2 + Y 2 



(40) 



The above complicated looking expression can be used to give a equally 
complicated expression for the gauge invariant momentum, however when 
plotted one obtains a very simple graph for the above. 

The graph in Figure (11) reveals that the behavior of the momentum is 
damped, but not drastically as the P r , also the singularity does not appear 
as any wild point here. This reaffirms the fact that the Momentum, being 
densitised functions of the dreibein should not be singular at the origin of 
the black hole. 

Thus we now have all the ingredients necessary to write the Coherent 
state wave-function for the planar graph. 



(39) 
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Figure 11: P t 
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Fine, so the momentum classically can be evaluated, and encodes infor- 
mation about the geodesic deviations which start dominating once inside the 
horizon, and inability of r = const surfaces to exist. But what about the 
coherent state, which is peaked about the above classical quantities? How do 
the fluctuations around the above classical quantities depend on the distance 
from the horizon? 



4 Quantum quantities 

The coherent state as defined in (5) here, can be written in a series in the 
irreducible representations. Using the parameterisation for h = e^ T \ where 
/j, 1 denote the angles or three SU(2) parameters, one can write down the 
series explicitly [17]. As per (5), this involves determining the irreducible 
representations and the traces corresponding to each for the group element 
gh^ 1 , and then finally summing over all possible representations. However, 
to see whether the function converges as t — > 0, i.e. when the width of the 
state is reduced, one has to use a Poisson re-summation formula. This is 
what was used in [17], and we quote directly from there. This form, given 
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the complexified element is given as: 



= t 3 / 2 fx^ t ^ ° (x)-2mn)exp^ 



2t 



n=— oo 



(41) 



Where 



x = cosh — cos u + % — — sinn — sinu 
\2 J P P9 V 2 J 

= cosh cos A 4 + * cos ( r ) sinh ^"2") s * n ^ 

= cosh(s + i(f>) (43) 

a; is the trace of the group element or x = Tr(g/i _1 ), and one can choose to 
write it in a compact form in terms of new variables s and <fi as cosh(s + i<fi). 

Using the above, one can in principle determine the expectation values 
of the various operators, defined on this generic edge. But first, let us study 
the peakedness properties of this state. Quoting [17], again we find that at 
t — > 0, the so called semi-classical limit, the probability goes as 

PUh) = ( 16 ;f l^'^V ^W) ( (44) 



\x 2 



with 5 2 = P 2 /4 — s 2 + 2 — fi 2 . Now this is the function of interest in 
the above probability amplitude. Clearly, as observed in [17], this function 
captures information about the non-abelian nature of the probability. Thus 
we concentrate on 5 2 now, as proved in [17], this function 5 2 > 0. Which 
also implies that as a function of 5 2 , the probability amplitude is maximum 
when 5 2 = 0. Where is 5 2 = 0?? Amazingly looking at the definition, this is 
zero when s = P/2 and (f> — //, which when looking at (43) is solved at 

cosr = ±l! (45) 

Let us thus look at how the above behaves as a function of the 'graph' 
degrees of freedom. The behavior of cos r thus determines the size of the peak, 
and also the sharpness of fall around the peak. Let us investigate along the 
following lines. Due to the fact that we have taken a Radial edge, and the 
particular gauge choice of diagonal connections, the classical holonomy is of 
the form h r = exp(/(r)T 1 ). Thus when we take the variable h = exp(//T 7 ), 
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then the /z 1 captures corrections to f(r), while fi 2 and /i 3 measure corrections 
along the T 1 and T 2 directions, where classically the holonomy component 
has no contributions. Now, let us examine (45), before moving on to other 
quantities: 

cosr = — — (46) 
X^r)^ + -y=5= [sin (7a 3 ) (a/12 + Ate) + cos (7a 3 ) (fi 2 - «Ate)] 

V x ! + x i yM + t4 + t4 

First the simple facts: When X 3 = 0, X 1 7^ 0, cosr = 1 when /j, 1 ' 2 = 0. 
Thus, when X 3 = 0, the most favored configuration is abelian, i.e. there are 
no corrections to the classical holonomy along the other two directions even 
at the quantum stage. However, the existence of such points in the classical 
phase space are not independent of the graph degrees of freedom, e.g. from 
(29,30), X 3 = whenever 9 = 7r/2, or the radial edge intersects the dual 
surface at a specified angle! However, this apparently trivial point acquires 
some meaning when one realizes that one is measuring distances still given 
by the classical metric, and, it is along 9 = n/2, one can think that the radial 
edge acquires the status of a line in flat space: sin 9 — 1, as the 2-surface on 
which the momentum is evaluated has metric d9 2 + d<f> 2 . In so me sense this is 
justified as one can think of this as the 'ray of least resistance' and as shown 
in (32), in flat space, the momentum is abelian for the radial edge. Figure 
(12) are the plots for P}jP r and P 2 / 'P r , while varying the angle at which 
the radial edge is positioned (6> = d) and the position of the radial edge 
itself embodied in a. Now there can be more solutions to cosr = 1, and each 
will come with it's own physical interpretation in terms of the classical black 
hole metric. Thus, the coherent state is very sensitive to classical geometry, 
and this is precisely what we were seeking. 

Both the plots in Figure(12) have been evaluated at 9 = 7r/100, i.e. with 
a very fine graph indeed, and confirm our claims with all the points at which 
Xi = 0,X 3 7^ and X x 7^ 0,X 3 = 0, appearing at 9 = 0, n and 9 = n/2 
respectively. However note that in the first case one does not get 5 2 = 0. 

Next, we show the plots for as one makes the dual surfaces smaller and 
smaller with the radial edge stationed at 9 = n/6 Figure(13). Here, P} / P r 
is plotted as a function of the size of the dual surface and also the position 
at which the dual surface is intersected by the edge, measured as a function 
of a. Oscillations set in inside the horizon, indicating that the solution is 



27 




Figure 12: P}.jP r P?/P r 

no- longer stationary, (r^/aw = 200). Note that in all the discussions above, 
the amplitude is moded out of the P giving finite values for all a. Thus 
the dependence in cos r on P r is not abnormally damped, and one can study 
interesting behavior for all values of a. The wavefunction is now completely 
determined, and one has to evaluate the expectation values of operators in 
these states, as a function of the distance from the horizon. 
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Figure 13: P}jP r 

This is a similar plot for P r 2 /P r , the interesting aspect of this plot was 
that without the modulating factor sin(7a 3 ), Mathematica gave a wild graph 
at a = 1, which prompted us to think, that we had discovered some new 
quantum gravity. The graph is attached below as an example of the many 
false alarms. 
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Adding the proper modulating factor for the P r 2 /P r , P r 3 / ' P r we as expected 
get good behavior Figure (14). (Note for some reason in Fig. 13 and 14, 
Mathematica generated 6 as q in the .eps file, and the P r s appear as 
should be read as p( 1,2 ' 3 ).) 
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Figure 14: P r 2 /P r and P?/P r 

All these graphs above, show the existence of certain points at which , 
the functions become highly oscillatory, and these points occur always inside 
the horizon with their distance decreasing from the centre, the finer and finer 
the graph is. Note that making the graph finer allows us to probe shorter 
and shorter length scales. Also note that since one has moded the amplitude 
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out, the functions have values which are quite reasonable. 

Ok, so much so good. But we wanted to calculate corrections to geometry, 
since anyway we know that the coherent state has the peakedness property 
as defined above. To do so, we take again the path of least resistance, i.e. 
take the radial edge for whom the coherent state behaves as a Gaussian. 

Figure (15) is a plot for cosr = 1 and 5 2 = 0, with the peak appearing as 
one takes t — > 0, precisely at /i = 
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Figure 15: Probability Amplitude: Peak 

Note the probability amplitude is not the probability and hence not re- 
stricted below 1. Looking at the probability amplitude, it has the following 
form as a function of P and fi. 

{ > W 3/2 (sinh 2 (P/2) + sin 2 Gu)) J\ P ) 1 ' 

This is almost gaussian in fi, except for the modulating factor in front. To 
examine that, we plot the above as a function of (P, /i) Figure(16). and 
t = IQ(~ 6 \ The plots show very good behavior. Next, let us look at /i = or 
the location of the peak in the probability amplitude. One finds this strange 
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behavior that the peak takes the P dependent value [18]. 

Psinh(P) 

(cOsh(P) - l)t 3 / 2 1 J 

This is rather surprising, which means that the wave function is peaked higher 
and higher as P increases, which is fine for the case of the black hole, where 
P increases polynomially in the radial distance outside the black hole. The 
increase in the probability amplitude implies that the wavefunction is more 
'semiclassical' in some sense. 




Figure 16: Peak Outside Horizon 

How about the behavior near the point where P oscillates. This is as 
illustrated in Figure (17). 

The first plot is for the P, and the next one for the wavefunction around 
that value of P. Note that this is a plot for black holes with size greater than 
10 9 /p. Hence as we discuss in section 4.2, the momentum can be assumed as 
continuous. 

So, apparently, everything is fine inside the black hole, along this radial 
edge. The coherent states appear very robust indeed. 
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Figure 17: Momentum Inside and Corresponding Peak 

This now brings us to the calculation of expectation values of the op- 
erators. Do the quantum fluctuations behave nicely? Or there are certain 
regions in the black hole where they grow to be very large? The corrections 
ofcourse shall be proportional to t, as computed in [17], but they might be- 
come uncontrollable in some region, signifying the breakdown of semiclassical 
approximations. To determine the corrections, we see the evaluation of the 
expectation values of the momentum operator as given in [18]. The expec- 
tation values along the edge of least resistance will have corrections of the 
form: 
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y 2 6 2 T(P 2 -P/2y + y 2 ) (cosh P cosh y/2 - sinh P sinh y/2)~ 
6XP ^4t ~ P [ (sinh P cosh y/2- sinh y /2 cosh P) 2 _ 



/ sinh P(sinh P cosh y — cosh P sinh y) , . 

x i / i — y —— ^ 49 

V p(p-v) 1 ; 

Where we have defined y = P' — P. This is done to extract all the P 
dependence out such that one is left with the integral in y, which being a 
dummy variable shall not yield any new dependence on P. This gives for 
large P corrections 

AP = t(P + 1 + 1/P) (50) 

For small P 



AP = t(P 2 + P + l) (51) 

Thus at least to order t, the fluctuations are perfectly under control. Can 
there be non trivial results for higher order corrections in tl At present we 
cannot have an answer on it. However, it is quite a miracle that P is every- 
where well behaved, otherwise, existence of such regions would mean break- 
down of the semiclassical approximations. The corrections to the holonomy 
operator as determined in [18] go as t times a polynomial of the momenta, 
and hence are finite everywhere. Any operator, be it curvature, are or vol- 
ume can be expressed in terms of the elementary operators, the holonomy 
and the momenta. Hence, quantum fluctuations for them will also be under 
control. 



4.1 A Quantisation Scheme 

It is interesting to observe that to study the full black hole slice, one need to 
use the Electric field representation. The graph in the configuration space had 
to exclude the singular point, but the dual graph by construction includes r = 
0. Also, as per equation (33), the singularity at the origin of the black hole is 
not seen by the dual momentum as some divergent region. Hence, to probe 
this region, one uses the Electric field representation. The wavefunction in 
the Electric field representation has been obtained in [17]. Following [19], 
the wavefunction in the 'electric representation' corresponds to coefficient 
of one representation 7ij imri of the Peter- Weyl decomposition. Instead of a 
continuous valued variable /i as in the configuration space, one obtains the 
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wavefunction as a function of j, mn corresponding to the Eigenvalues of the 
momentum (j corresponds to the eigenvalue of the Casimir P, while m and 
n correspond to the eigenvalues of the left and right invariant P3 [17]). The 
wavefunction is simply: 

^(jmn) = e-W+^ig)^ (52) 
The probability appears as: 

p. Kexp {_Ki±v^!!} (53) 

Where P is the classical gauge invariant momentum, and j gives the 
eigenvalue of the Casimir in that particular representation. Clearly for jt > 1 
one can obtain a continuous value for P. But as mentioned earlier, this 
is not true when the momentum becomes very small indeed. Let us take 
an example, as evident from previous section the semiclassicality parameter 
being lO* -9 ) is enough to obtain the peak. Whereas inside the horizon, the 
Momentum oscillates reaching these values very quickly [Figure(8)]. The 
expectation value of P from Equation (53) is (j + l/2)t, and for one radial 
edge, one can determine a unique j. However, it also means that if there is 
an adjacent radial edge, with a momentum value close to it, then P r — P' r = 
(ji — J2)t. The discrete nature of the spectrum P = (j + l/2)t emerges and 
only those classical values are 'allowed' as peaks in the Coherent State. 

No-matter how big the black hole is, once the oscillations start, the ampli- 
tudes decrease very quickly, and soon, one is measuring values of momenta, 
which are 0(t) even at their zeroeth value, and hence here, comes the inter- 
esting conclusion, that the Peak values must be discrete!. 
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(54) 



Most notable is the presence of the remnant value 1/2 which prevents 
the momentum value from ever being zero. This shows that the point r = 
where the momentum necessarily vanishes, must be excluded from the 
quantum spectrum. Thus if the black hole is in a coherent state as defined 
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in [18], then only discrete values of the classical momentum can be sampled 
by the expectation values. What can we conclude? Quantum mechanics 
dictates that near the singularity, radial edges can be embedded only at 
specific quantised values. Does that mean that space is actually discrete? 
Will the above quantisation lead to an entropy counting? It is too early to 
say all that. However it is a certainty that the singularity will be smoothened 
out, due to the quantum uncertainty principle. 

5 The Actual Graph 

Ok, so much so for the single spoke, but it was needed to set direction to 
the calculations, and what quantities to look for, when one considers the 
actual planar graph. Let us begin by adding more equally spaced spokes, 
and corresponding dual edges. The actual graph will have tiny radial edges, 
tiny spherical edges along 9 and <p directions, and the dual surfaces are cor- 
responding two surfaces, stretched transversely to the edges. Now, given a 
vertex of the graph, there will be three edges starting out along the radial, 
and the two angular directions and similarly three going in. However, as an 
act of simplicity, the angle will again be suppressed. So the coherent state 
at each vertex will be 

— ^Pe i r n 4'e i g n 4'e° ut 4'e^ ut (55) 

As before, the non-trivial information must be in the ip er , and we proceed to 
extract that. The radial edges are now joined end to end, with dual surfaces 
intersecting them at finite distances from each other. 

To determine the momenta (8), one has to look at the surfaces, which 
make up the sphere now. Choosing a regular lattice/graph, one gets these 
little dual surfaces of equal length for a given radius. The radial edge inter- 
sects them at the centre transversely. Thus on a plane, the cycle consists of 
N vertices joined by edges of length e/2 and one unit of dual edge comprising 
of two such little edges. The vertex in the midpoint is actually the point of 
intersection of the radial edge with the dual edge. This labelling helps in 
evaluation of the holonomies and the momenta for the graph. As illustrated 
in Figure (18), one starts measuring the distances along the surface from 
= 0, and the surface is divided into N bits of length re/2 each, such that 
Nre/2 = 2nr. The bits are numbered starting from 9 = point as m = 0, 
and then + e/2asm=l and so on. The radial edge intersects the little 
surface at m = 2n + 1 odd vertices (denoted in blue dots in the diagram), 
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and the little surfaces extend from m = 2n to m = 2n + 2 vertices(denoted 
by red dots in the diagram). The value of the radial holonomy is as given in 
the previous section Equation (24). For the calculation of the momentum, 
one needs the holonomies on the surfaces, which are given by (25), where 9q 
is the point where the radial edge intersects the two surface, and is denoted 
by 9 = (2n + l)e. 




m +2 



m+1 



Figure 18: The Labelled Vertices 



The momentum are now given by 



1 



P l e = Tr [7* he {X{r))h~ l ] 

UN 



(56) 



With X(r) being evaluated for the range of 9 = 2me to 9 = (2m + 2)e. which 
gives 



Xf(r) 



r 

T 



sin(2m + l)e 



sinfl — a')e sin(l + a')e 



(!-«') 



(1 + a') 



(57) 
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Since we have already dealt with the behavior of the above, we donot go 
into detailed analyses again. The only thing to take care about is that the 
radial edges are stationed at discrete values of 9, and the size of the dual 
surface is fixed at 2n/N, N' denotes the number of the radial edge (N' will 
be same for radial edges at same distance from the center) . Thus the only 
variable is now N given N', and we can choose it to be as large as possible. 
The momenta will be oscillatory at greater values of a as N is increased. 
Hence finer and finer graphs probe shorter and shorter length scales. 

Also for the angular edge, we again determine P/, as determined earlier in 
Equation (38,39), however now the limits of integrations are from (2n)S/2 to 
(2n + 2)<5/2, with the angular edges intersecting the dual edges at (2n + l)S/2 
radial distances. The dual graph obviously has the zeroeth edge which starts 
out from the origin. Though the fact that the individual components in 
(38,39) look singular at the origin, the total expression for the gauge invariant 
momentum P = \/Y^ + F 2 2 is just given in terms of a + and a_. Taking a_ = 
and a + = 7, one obtains finite terms in the expression given in Appendix 
3. (Lim^oo sin(x)/x = 0). Thus, one can write down the Coherent state for 
the entire planar graph. 

Since the individual coherent states are peaked at their respective values, 
the entire graph will also be so. The non-trivial expectation values has to 
come from the evaluation of Volume or Area operators in the Black Hole 
Back ground. This calls for a entire new calculation [48]. Here, we make the 
following observations. Suppose now, we look at a single radius r, and take a 
product of the coherent states of all the spokes of the wheel which intersect 
the radial surface. This should be a 'spherically' symmetric coherent state. 
Thus, it would be: 

N 

As long as one is not evaluating expectation values of operators linking 
various edges, one can take the final probability amplitude to be a product 
of the individual probability amplitudes: 

p\r)= \[p\h n ) (60) 

N 

And hence one can define an effective probability amplitude, with an 
effective momentum etc. The effective momentum goes as l/a AN , depending 



cos(2m + l)e 



sin(l — a')e sin(l + a')e 
(1 - a') + 1 + a' 



(58) 
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on the number of radial edges one introduces. Hence, as the graph is made 
finer and finer, a > 1 regions, the value of the momentum decreases rather 
fast. The oscillations will now be with a effective size of the sphere which is 
almost macroscopic and set in very near the horizon. In some sense one will 
be soon in a region of microscopic physics for small black holes. The details 
are under investigation. 

6 Discussion 

In this article, we addressed the coherent state for a black hole. Our main 
contribution is the evaluation of the graph degrees of freedom h e , P e for a 
spherically symmetric graph in a constant time slice of the Schwarzschild 
black hole. The coherent state by construction is peaked at the classical 
values of these, and we observed that at certain radial edges, coinciding with 
radial geodesies, one gets almost Gaussian behavior of the wavefunction. In- 
side the horizon, the geodesic deviation forces become very large, and the 
r = const spherical surfaces are torn apart, which reflects in oscillatory be- 
havior of the momentum P/. These can be thought of as regions of failure 
of the graph embedding, and as one uses finer and finer graphs these regions 
recede into the black hole. Do we get a smoothening out of the black hole sin- 
gularity? The quantisation condition on the momenta for small black holes, 
suggest that indeed space is discrete inside the black hole, in the vicinity 
of the singularity. How robust is the result? How do the graph degrees of 
freedom affect these, when one tries to average over them? One does not 
know, but surely, a very strong indication emerges of the physics to come. 
As of now, we discuss the following points: 

• The only other 'non-perturbative' semi-classical treatment which in- 
volves calculations to black hole geometry exists in 2-dimensional CGHS 
model, though it has matter coupled to it [29] . Perhaps, a 2-dimensional 
version of these coherent states can be written down, and the results 
compared with existing literature. A canonical framework for the 2- 
dimensional model already exists [30]. 

• The quantisation conditions which arise seem to suggest a degeneracy 
associated with the black hole, and one can get a calculation of Black 
hole entropy associated with the entire space-time through the coherent 
state. Will it in that case, reveal what the Hawking radiation spectrum 
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should be [31]? All these and many more issues involving corrections to 
entropy, [32, 33] are yet to be addressed but the framework now exists. 

The value of (3 gets fixed here if one comments on the classical hori- 
zon at r g = r, however taking into account the quantum corrected 
geometry, one may find (3 to be the same value predicted from pre- 
vious approaches, including the recent quasi-normal mode approach 
[34, 35, 36]. The details are under investigation [48]. 

Classically, the holonomy and the Momentum variable have oscillatory 
behavior, as one approaches the singularity. These are non-local ob- 
jects, and donot correspond to the usual oscillatory approach to singu- 
larity associated with scale factors in cosmological models. In ordinary 
Schwarzschild black hole, the approach to singularity is monotonic [37] , 
however slight perturbations lead to a Oscillatory approach. There is 
a possibility, that the Momentum e.g. can be written down as f *E', 
where now the E corresponds to a different metric, (£" has all the 
holonomies associated with the Schwarzschild absorbed in it), in which 
case, the oscillatory approach can be associated with a BKL treatment 
of singularities [38]. 

The original derivation of Hawking radiation involves the Event hori- 
zon, which is globally defined. However, here, the knowledge of few 
adjacent slices and the presence of the apparent horizon is the only 
ingredient. Unless evolution under the Hamiltonian is ensured, one 
should be careful. BF theory on a single slice has the same results as 
Einstein Gravity [39] , and hence it is important to ensure that inclusion 
of time evolution does not change the results drastically. 

There exists semiclassical derivation of the expectation value of the 
area and volume operator [40, 41, 42] in [46, 47]. The same techniques 
can be used for the black hole, and corrections to the area operator 
should be particularly interesting. 

One must couple matter to gravity and repeat the same exercise as 
most semiclassical derivations of Hawking radiation and back reaction 
involve matter fields, and their super-planckian energies in the vicinity 
of the horizon [43, 44, 45]. Matter coupled semiclassical states are 
treated in [46, 47], and they need to be extended for quantum fields in 
the curved space-time of the black hole. 
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Though the results reported here need more concrete interpretations, as a 
begining, there has been many interesting observations, and many avenues 
revealed which need to be explored. 
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Appendix 1 

The actual graph consists of a foliation of 3-dimensional spheres, and their 
appropriate duals. The dual surfaces are 2-dimensional, and involve both the 
4> and 9 integrals. The radial edge thus intersects a 2-surface which is part of 
a sphere sandwiched between two such spherical embeddings of the original 
graph. This graph is effectively three dimensional and inclusion of the <fi 
coordinate adds some complications, but does not add anything drastically 
new to the observations already made in the previous section. As usual 
we begin with a single radial edge, and calculate the corresponding three 
momentum now. The 'dual polyhedronal decomposition' of the spherical 
geometry now involves pieces of spheres, eg, the dual surface to the radial 
edge is piece of a 9, 2-surface. This ensures that the radial edge intersects 
it transversely and only once. The projections of the dual surface onto the 
plane, gives the edges of the dual planar graph discussed in the previous 
section. There is however a subtle ambiguity associated in choosing the 
edges of the dual surface which are used to convolute the triad. One way to 
fix this is to embed the edges so that they lie along 9 and directions, and 
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then define the path from any generic point to the point of intersection to be 
along the shortest distance curve. However as discussed in the Appendix 2, 
this method leads to singular answers for the momentum, and is not suitable 
in the study of this problem. One thus resorts to the choice of a path from a 
generic point 9, to the point of intersection Q , O to be along a = path 
which takes 9 — > 9 and then moving along 9 = 9 curve such that — > O . 
The momentum evaluated thus should have now the following integrations: 

P e f m = - — Tr[TX / V* * Eh?h?h£\ (A-l) 

UN J 

Surprisingly X 2 = 0, even in three dimensions, and the origin of that lies 
in the choice of the point of intersection to be exactly in the middle of the 
dual slice. The X[ and the X' 3 are related to the 2-dimensional ones by the 
following simple relations: 

X[ = eXi (A-2) 
X' 3 = eX 3 (A-3) 

e is the width in the direction, and it is clear that for e small the above re- 
duces to (29,30). The invariant momentum is different from the 2-dimensional 
one in it's dependence on sm(neasin9 ) , and has very similar properties. The 
interesting aspect is that the momentum never depends on the position of 
the radial edge in the direction, which is a manifestation of the azimuthal 
symmetry. The momentum for the angular edges is modified in a similar 
manner, and one obtaines the coherent state, as a product of the coherent 
states for all edges. 

Appendix 2 

In this section, we discuss the ambiguities associated in the choice of edges 
on the 2-surface dual to the radial edge. The most appropriate choice of 
edges would involve the choice of the shortest path connecting a generic 
point 9,(j) to the point #O)0o- This, would involve if 9 < ir/2, and 9 < 9 , 
a path along 9 = 9 such that — ■> O an d then a path along = O and 
9 — > 9 . The situation will be reversed for points 9 > 9 , and one has 
to travel along = path to 9 — > 9 and then along 9 = 9 . This en- 
sures that the contribution to the invariant distance comes from min sin 9. 
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Thus, one has to evaluate the integrals of the 2-surface points 0, sepa- 
rately for < 9 and > 9 . Thus for 9 > O , h(9,<f>) = h{<t>)h{9) 



cos 



^0 I 



i - 0)) - ta 3 sin(f (0o - 0))] [cos(f(0 o - #)) + ^ 2 sin(f (0 O - 0))] and 



for < O , this is h(0, 0) = h(0)h(<j>) = [cos(f(0 o - 9)) ~ ™ 2 sin(f (0 O - 0))] 
[cos ( | (0o — 0)) — «cr 3 sin(|(0 o — 0))], (z/ = a sin O ). The Integrals (Note for 
simplicity all calculations are done for a 3> 1: 



h(0,0) *Eh(9,(f))- 1 
= r 2 J \fo-ia I f I )a 1 {f + ia 1 f T ) sin, 
= [(l- 2/ 2 2 - 2/ 3 2 ) ( x 1 + (/o/3 + /2/i)^-(/o/2-/i/3)^]sin^0 

= r 2 y (cos a(0 o — 9) cos z/(0 o — 0)<r 1 + sin z/(0 o — 0)<x 2 — sin a(0 o — 9) cos z/(0 o — 0)o" 3 ) sin 9d9d> 
= r 2 J (cos a(9 — 9) cos ^o(0o _ 0) " 1 + sin ^o(0o ~ 0) c °s cy(9 — 9)a 2 — sin a(9 — 0)cr 3 ) sin 9d9i 

The interesting observation is the fact that the integral proportional to a 2 is 
zero for both. The integral proportional to a 1 yields the following for 9 < 9 
(one takes as the limit: O — ne/2 < < o +ne/2, O — me/2 < 9 < 9 +me/2): 



r 
a 



^ J 2 fe+i(o;ne) I cos[(2/c + 1)0 O 



fc=0 



2k + 1 



(2fc + l) 2 -a 2 



(A-5) 



cos[(2fc + 1)0 O — (2Jfe + 1 — a)me] cos[(2A; + 1)0 O - (2k + 1 + a)me] 



2A; + l-a 2A; + 1 + a 

Thus the combined Pi component now has the form: 



2r 2 

Pi = ' 



cr 



(2k+l) 2 -a 2 



J ^k+i{ane) [{2k + 1) cos(2A; + 1)9 



(A-6) 



k=i 



x 



(2k + 1) cos[(2A; + l)(0o — me)] cos(me) a sin[(2A; + l)(0o — me)} sin(o;me)}] + = 

1 or 

sin(naesin0 o )\ a n , , , ^ ■ a ■ ( \ ■ ( \ 

Jo(ane) cos0oU _ cos(me) cos(amej) — asin0 o sm me sin fame) 

sin O / 



, . sin(anesm9 , , . . 

+ Jo(a>ne) H ( sm O sm( /le) cos( «me + a; cos O sm(me) smfamej) 

sin O ' 
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However complicated the above expression is, the most important observa- 
tion is that the expression for momentum is singular at a — 2k + 1 , and this 
is not a good sign. The moment one chooses the paths which have symmetric 
behavior around the point of intersection, one gets smooth answers. In oth- 
erwords, if one chose to consider the choice of path for 9 < 9 and extend it 
to 9 > 9 without thinking about the shortest possible way, then one would 
get for Pi the following: 



p _ 9 



^ J 2 k+i(^ne) sin(2A; + 1)9 



k=0 



sin(2/c + 1 — a)me sin(2A; + 1 + a)me 



2k + l-a 2k + l + a 

(A-7) 

Note in the above, one gets higher modes of oscillations, but the nature is 
same, the oscillations set in after a = 1 and are highly damped inside the 
horizon. Thus, the details of the graph, and edge choices drastically change 
the nature of the oscillations, what then is the graph invariant physics? Obvi- 
ously, one needs to be careful about the embedding in the above space-time. 
However the relief is that, making the graph finer and finer, such that argu- 
ments are smaller (ne) and (me) the ambiguities are reduced even for large 
values of a, such that one keeps only the Oth term in the above series. 



Appendix 3 

Here we quote the gauge invariant momentum for the angular direction. We 
have relegated it to the appendix due to the length of the expression. 



Pee = 



t 4 sin 9o 

a,N 



+ 6a l) 2 + (_ Q 2 + 6 « 2 ) 2 + {-af + 2af ) ' 



+ {-ai 2 + 2aJ 2 y + [Ctia 1 / 2 ) - Ct(aJ 2 )j ' + [St(a/ 2 ) - Si{a! 2 ) 



2 cos 



1/2 1/2 ) V 

a_ a+ J K 



-a + + 6a 2 h )(— a_ + 6a 2 ) + (- 



-a 



1/2 



2a 3 / 2 ){-a 1 ! z + 2a 



1/2 



+ 2sin(4^ - ~Te) {(-«+ + 6a 2 h )(-a V2 + 2a 3 / 2 ) - (-a_ + 6a 2 _)(-a_/ 2 + 



2a 



3/2 



- 2[Ci(a l l 2 ) - Ci(at! z )} <j cos(^-)(-a + + 6a 



1/2 



a 



V2' 
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cos 



a 



1/2 



)(-«_ +6a 2 ) 
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+ (-a 1 / 2 + 2af ) sin(-L) - (~«- /2 + 2a 3 _ /2 ) sin(-^) 

a + 7 a_ 



2[^(af ) - Si(a]/ 2 )} \ cos(4^)(-«l /2 + 2af) + cos(-^)(-ai /2 + 



1/2 ^ 2a S/2x 



+ sin(— j-^)(-a + + Qa\) - sin(— j-^)(-a_ + 6a 2 ) 



1/2 
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